Period tripling causes rotating spirals in agitated wet granular layers 
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Pattern formation of a thin layer of vertically agitated wet granular matter is investigated exper- 
imentally. Rotating spirals with three arms, which correspond to the kinks between regions with 
different colliding phases, are the dominating pattern. This preferred number of arms corresponds 
to period tripling of the agitated granular layer, unlike predominantly subharmonic Faraday crispa- 
tions in dry granular matter. The chirality of the spatiotemporal pattern corresponds to the rotation 
direction of the spirals. 

PACS numbers: 45.70.Qj, 45.70.-n, 45.70.Mg 



From the formation of galaxies to hurricanes, from the 
structure of a seashell to the magnetic ordering on the 
atomic scale pL], spirals are ubiquitous in nature [2]. Dur- 
ing the past decades, various spiral forming systems have 
been extensively investigated, e.g. Rayleigh-Benard con- 
vection Q, the^ Belousov-Zhabotinsky reaction [4], and 
cardiac tissue [5]. Moreover, spacecraft observations of 
the Saturn's rings - granular matter consisting of icy 
particles - reveal that the perturbations arising from the 
gravitational field of its moons propagate in the form of 
spiral waves [6], which have been modeled in terms of a 
hydrodynamical description 0. Despite this success, a 
continuum description of both dense and dilute granular 
systems [8] , including wet granular systems [9] is still far 
from complete. 

In this Letter, we present rotating spiral kink waves 
formed in agitated wet granular layers. We demonstrate 
that the dominating three armed spiral is a manifesta- 
tion of the period tripling in the system, and that the 
spatiotemporal chirality associated with period tripling 
drives the rotation of spiral arms. In contrast to those 
in dry granular systems [10], the arms are kinks sepa- 
rating regions colliding with the container at different 
phases. Understanding this particular spiral forming sys- 
tem could not only provide a step towards a continuum 
description of dry and wet granular matter, but also 
broaden our knowledge about the more general spiral 
forming systems mentioned above. 

Fig. [1] shows a sketch of the experimental setup. The 
granular sample is prepared by adding purified water 
(Laborstar TWF-DI) into glass spheres (SiLiBeads S) 
with a diameter d = 0.78 mm and 10% polydispersity. 
They were cleaned subsequently with ethanol, acetone 
and purified water, and dried in an oven before use. The 
liquid content W = V^/Vg, where is the volume of the 
water and Vg that of the glass beads, is kept within a few 
percent. A cylindrical polycarbonate container with an 
inner radius R = 8 cm and a height 1.05 cm is mounted on 
the electromagnetic shaker (Tira TV50350). The amount 
of sample (mass m) added corresponds to a height of a 
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FIG. 1. (color online) The experimental setup with sample 
images of a spiral pattern. Vertical or horizontal information 
of the patterns formed is obtained by a high speed CCD cam- 
era mounted above the sample illuminated by a laser sheet or 
a low angle LED ring correspondingly. The camera is trig- 
gered by a multi-pulse generator. Snapshot (a) is captured 
with the laser sheet illumination. The red half transparent 
line is a surface profile determined by an image processing 
procedure. The spiral arm corresponds to a kink between 
regions with different heights. Snapshot (b) corresponds to 
the low angle illumination (averaged over three synchronized 
images to enhance the contrast). 



few particle diameters. The sinusoidal shaking frequency 
/ and amplitude of the shaker are controlled by the func- 
tion generator (Agilent FG33220). The nondimensional 
acceleration F = An'^pA/g^ with vibration amplitude A 
and gravitational acceleration is obtained by an ac- 
celerometer (Dyson 3035B2). 

A high speed camera (IDT MotionScope M3) mounted 
above the container is used to capture the patterns: the 
LED ring illumination provides the horizontal motion of 
the patterns, while the surface profile is obtained via the 
laser profilometry method [11]. The laser sheet illumi- 
nates a line on the sample surface with an angle = 38°. 
The camera is triggered by a synchronized multi-pulse 
generator to capture images at fixed phases of each vi- 
bration cycle. 
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FIG. 2. Snapshots of migrating fronts (left panel) and a ro- 
tating spiral (right panel) captured at a fixed phase of con- 
tinuous vibration cycles (indicated as numbers) showing the 
period tripling behavior. Parameters: / = 80 Hz, W = 1.6%, 
m = 96g and r = 36(left), 37(right) 



Typical patterns observed are rotating spirals, which 
coexist with migrating fronts occasionally generated from 
the side wall or nucleated within the bulk (see Fig. [2]). 
One striking feature arising from the snapshots is that 
both patterns show a period tripling behavior: similar 
patterns appear at every third vibration cycle, in con- 
trast to the subharmonic patterns commonly seen in agi- 
tated dry granular matter [12., J^] . The migrating fronts 
tend to merge with one of the spiral arms and vanish, 
because both are of the same origin: they are kinks sepa- 
rating regions of different phase shift within the period-3 
vibration, as will be demonstrated below. The rotation 
period for the spiral shown in the right panel is ~ 1.8 s, 
much larger than that for the driving 1/80 s. 

There is a clear threshold for the patterns to emerge, 
as the step change of the order parameter I^^ns with F in 
the inset of Fig. [3] indicates, /rms = VYl {H^^ v) ~ 
with /(x, y) the image intensity at the position (x, y) and 
I the spatial average of /(x, y)^ is chosen as the order pa- 
rameter. No clear hysteresis could be detected within 
the error of F, as step increasing (Fine) oi" decreasing 
(Fdec) the acceleration for various waiting times At yields 
the same threshold, which is determined by the averaged 
steepest slope of /rms- Within the range of frequencies 
investigated, the threshold is increasing monotonically 
with frequency and insensitive to the geometry of the 
container and the liquid content W: measurements for 
/ = 80Hz with a square shaped container (with an in- 
ner side length 10 cm) and with a varying W up to 9.4% 
yield qualitatively the same spiral patterns and quanti- 
tatively the same threshold. The weak dependence on 
the hquid content, even for the regime where large liquid 
clusters form, can be attributed to the insensitivity of the 
rigidity of the wet granular material on W [14], which is 
important for the kink dynamics because the granular 
layer will be bent and stretched. Nevertheless, wetting is 
a necessary condition for the present instability because 
no patterns could be observed without wetting liquid. 
Below 55 Hz, visualization of the pattern is prohibited 
by the particles sticking to the lid. 

In the pattern forming regime, the number of rotat- 
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FIG. 3. (color online) Phase diagram for the patterns (shaded 
with squares) measured by varying F while keeping the vi- 
bration frequency for m — 129 g and W — 1.6%. The region 
with F > 50 is unexplored. The snapshot (averaged over 
three synchronized images) shows two spirals rotating in op- 
posite direction captured at / = 70 Hz and F = 23. The 
red dash curves highlight the spiral arms and yellow arrows 
denote the rotation directions. The inset shows the order pa- 
rameter /rms vs. F at / = 80Hz for various time step At 
and both directions of F variations. 



ing spirals increases with F. The spirals all show pe- 
riod tripling, and may rotate in the same or opposite 
directions (see the snapshot in Fig. [3]). They are con- 
sidered to be stable, since no decay could be observed 
for at least 15 minutes. Even higher F finally lead to a 
homogeneous gas-like state, where the particles occupy 
the whole space of the container. Decreasing F shows 
the reversed scenario: from a homogeneous gas state to 
the spiral patterns, until a flat featureless state appears 
below the threshold. In contrast to the phase diagram 
investigated before [15], no granular 'gas bubbles' could 
be observed here, the reason for which is a subject of 
ongoing investigations. 

To investigate the spiral dynamics, the surface pro- 
files are measured with laser profilometry. Fig. 2] shows 
the oscillations of an illuminated line close to the core 
of a rotating spiral arm at different time scales. The 
long time behavior (Fig. |4l^a)) indicates the propagation 
of the kinks separating regions with different heights H. 
Observations with combined illuminations indicate that 
the kinks correspond to the arms of spirals. This ex- 
plains why the arms look brighter: curved surface of the 
kink regions reflects more light from the LED ring to the 
camera. The birth of two counter propagating kinks cor- 
responds to the crossing event of one spiral arm with the 
illuminated line. Since the spiral core is close to the illu- 
minated line, the constant propagation speed suggests a 
linear relation between the radial distance to the spiral 
core r and time t, i.e. an Archimedean spiral, assuming 
a constant angular velocity. Once the arms are formed, 
they are resistant to the disturbances along the propa- 
gating arms. 

At the time scale of vibration cycles. Fig. IH^b) reveals 
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FIG. 4. (color online) Time space plots of the surface profiles 
obtained with laser profilometry. The long time behavior (a) 
presents the traveling of the kinks, i.e. spiral arms. The short 
time behavior (b), corresponding to the region marked with 
dash line in (a), shows the periodicity of the height fluctu- 
ations. Other parameters: / = 80 Hz, F = 32, frame rate 
320 Hz, H—0 chosen as the lowest surface height. 



the period tripling behavior of the height fluctuations, 
which is in accordance with observations from Fig.[2j The 
fluctuation amplitude being very close to the distance 
between the surface of the sample at rest and the lid 
7mm) indicates that the lid plays an important role 
in developing the spirals. In an open container, no such 
patterns could be observed. 

Fig. [5] is meant to illustrate the rotation mechanism of 
the spiral arms. Fig.[5fa) shows the spatially resolved co- 
variance between two subsequent images, which is calcu- 
lated by EAx,A2y-o y^Ay)It^i{x^Ax, y^Ay), 
where t denotes time and d the particle size in pixels. 
The time step (2.5 ms) is short enough for the mobility 
of particles to be captured (bright colors corresponding 
to smaller covariance, i.e. higher mobility). As inter- 
preted in Fig. Efb) and (c), the higher mobility of par- 
ticles (granular temperature) in region I arises from the 
fact that those particles have just undergone a collision 
with the container bottom. Region II and III will un- 
dergo this event in the consequent two vibration cycles. 
The bright line in Fig. [5fa) separating region II and III 
is due to the rotation of the spiral arm. 

The qualitative idea to explain the rotation of the spi- 
rals is the assumption that regions with higher granular 
temperature will expand into 'cold' regions. Thus we ex- 
pect a maximum velocity for the spiral arm separating 
region I and II because region II is the 'coldest' due to 
the longest time interval having elapsed after the last 
collision with the bottom. Therefore the predominating 
motion of the spiral arms will be the rotation towards 
a region that is going to collide with the container in 
the next vibration cycle. Reversing the time sequence 
among region I, II and III - corresponding to spatiotem- 




FIG. 5. (color online) (a) The covariance between two sub- 
sequent images taken as the collision with the container oc- 
curs. The brightness corresponds to the mobility of particles. 
Sketches (b) and (c) indicate the spatiotemporal chirality of 
the pattern. The black line in (c) indicates the positions of 
the vibrating bottom, the dash lines correspond to the center 
of mass in regions I, II and HI. H denotes the height. 



poral pattern with different chirality - would reverse the 
rotation direction of the spiral. Thus we expect rotating 
spirals for periodic states with multiples of the vibration 
period which are larger than 2. The migration fronts 
shown in Fig. [2] are explained by the same mechanism: 
the front moves towards the region lagging behind with a 
phase shift of 120°, corresponding to a time lapse of one 
vibration cycle. 

To get a clue why period tripling is so dominant in 
our apparatus, we use a single particle bouncing model 
to estimate the periodicity, which is the period scaled by 
1/f. It treats the sample as one wet particle colliding in- 
elastically with the container. In addition to model used 
to describe agitated dry granular layers [13], it considers 
the capillary force acting on the granular layer [9] and 
collisions with the lid of the container. Initially the par- 
ticle moves together with the container bottom until the 
vibrating acceleration is large enough to overcome the 
capillary and gravitational forces acting on the particle. 
After the first detachment, the particle takes a parabolic 
flight until the next collision occurs. After that, it may 
detach immediately or stay with the container, depending 
on the colliding phase. This process will run iteratively 
until the particle rests on the container bottom again, 
which ends the cycle and determines the periodicity. Ex- 
perimentally the periodicity is determined by fitting a 
Gaussian profile around the peak of the surface height 
fluctuations in Fourier space. 

Fig. [6] indicates that the simplified model captures the 
dependance of periodicity on F and predicts the thresh- 
old for time dependent patterns. Periodic oscillations 
appear above a certain threshold, where the period is an 
integer multiple of the driving period. Increasing F leads 
to larger periods because the free flying time increases 
as long as the particles do not touch the lid of the con- 
tainer. Once they touch the lid, the periodicity starts 
to decay due to the additional impact from the lid. Dif- 
ferent from the model, experimental results show high 
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FIG. 6. (color online) Periodicity as a function of driving ac- 
celeration for / = 80 Hz, obtained experimentally from the 
surface height fluctuations (a), and estimated by a single par- 
ticle bouncing model (b). Errors in (a) correspond to the 
variance of the Gaussian profiles fitted to the spectrum of 
the height fluctuations. Red solid symbols highlight period 3 
regimes above the threshold. Other parameters are the same 
as in Fig|3l 



uncertainties of periodicity below F ^ 30, which arise 
from the intermittent fluctuations with small amplitude. 



For F > 30, both the experiment and the model reveal 
a wide range with period tripling, which corresponds to 
the time dependent patterns observed in the experiment. 
In this regime, the periodicity is robust against small 
surface fluctuations. Note that the wide period tripling 
region appears after impact with the lid, which can thus 
be considered as a technical trick to expand the existing 
range of spirals. 

In conclusion, three armed rotating spirals are found 
to be the dominant pattern in agitated wet granular mat- 
ter. The preferred period tripling could be understood 
by a simplified model considering complete inelastic colli- 
sions between the granular layer and the container. The 
spatiotemporal symmetry breaking connected with the 
period tripling explains the rotation of the spiral arms. 
A more detailed model will be needed in order to calcu- 
late the exact curvature of the spiral arms. In fact, this 
pattern could provide a convenient test of more elaborate 
theories for the dense flow of wet granular matter. 
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